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Abstract. A short survey of the renormalization problem in QCD and its non-perturbative solution by 
means of numerical simulations on the lattice is given. Most emphasis is on scale dependent renormalizations, 
which can be reliably addressed via a recursive finite-size scaling procedure employing a suitable intermediate 
renormalization scheme. To illustrate these concepts we discuss some — partly recent — computations of 
phenomenologically relevant quantities: the running QCD gauge coupling, renormalization group invariant 
quark masses and the renormalization of the static-light axial current. 



1 Introduction 

Apart from its well established role as a non- 
perturbative framework to calculate relations be- 
tween Standard Model parameters and experimental 
quantities from first principles jij , Lattice Field The- 
ory is particularly designed to solve various renormal- 
ization problems in QCD j2[ ||. Since renormalized 
perturbation theory as analytical tool is limited to 
high energy processes, where the QCD coupling is suf- 
ficiently small, but inadequate for bound states and 
momentum transfers of the order of typical hadronic 
scales, fi ~ lGeV/c, a genuinely non-perturbative 
solution of the theory is generally required. This is 
achieved by numerical Monte Carlo simulations of the 
Euclidean QCD path integral on a space-time lat- 
tice. Though renormalization is an ultraviolet phe- 
nomenon (relevant scales ~ a) and QCD asymp- 
totically free, tolerable simulation costs prevent the 
lattice spacing a from becoming much smaller than 
the extent of physical observables so that a truncation 
of the lattice perturbative series is often not justified. 
Therefore, it is far more safe to perform renormaliza- 
tions non-perturbatively. 

In addition, Lattice QCD has a large potential to 
address the computation of fundamental parameters 
of the theory, which escape a direct determination by 
experiments. The most prominent ones among them 
are the QCD coupling constant itself and the quark 
masses, whose running with the energy scale is desir- 
able to be understood on a quantitative level beyond 
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perturbation theory — the central subject of the next 
sections. The knowledge of these quantities (e.g. at 
some common reference point, a s (Mz) or m(2 GeV)) 
might then also provide essential input to theoretical 
analyses of observables of phenomenological interest. 
For instance, the mixing ratio e'/e in the neutral kaon 
system incorporates the strange quark mass value. 

The renormalization properties of many other 
quantities have been investigated with lattice meth- 
ods, e.g. (bilinear) quark composite operators, AS = 
2 matrix elements and, as presented at this confer- 
ence too, structure functions [Q. Later we will briefly 
discuss the non-perturbative renormalization of the 
static axial current as a further example. 

During the last few years the lattice community 
has seen much theoretical and numerical advances 
|l], Here it is worth to mention at least the is- 
sue of O(a) discretization effects inherent in the Wil- 
son fermion action. In case of the quenched approx- 
imation to QCD, where all dynamics due to virtual 
quark loops is ignored, they have been systematically 
eliminated through a non-perturbative realization of 
Symanzik's improvement programme J^, |^, [|. Hence, 
lattice artifacts can be extrapolated away linearly in 
a 2 , which allows to precisely extract many physical 
quantities in the continuum limit, a — > 0. 

2 Intermediate schemes 

As a representative example for a non-perturbative 
renormalization problem we may consider the calcu- 
lation of quark masses through the PCAC relation, 

-FkTOk = (m u + m s )(0|u7 5 s|K) (1) 
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( U T5 S )mS = ^p(50,dA t )('"75s)lattice, (2) 

in which the scale and scheme dependent renormal- 
ization constant Zp relates the lattice results to the 
MS scheme and is computable in lattice perturbation 
theory. But since this expansion introduces errors 
which are difficult to control, a non-perturbative de- 
termination of the renormalization factor is needed. 
A non-perturbative renormalization condition be- 
tween the two schemes can, however, not be formu- 
lated, because MS is only defined perturbatively. 

The idea to overcome this problem is the intro- 
duction of an intermediate renormalization scheme: 
the lattice observable is first matched at some fixed 
scale fiQ to the corresponding one in the intermedi- 
ate scheme, and afterwards it is evolved from (j,q up 
to high energies, where perturbation theory (PT) is 
expected to work well. Nonetheless, as in a simu- 
lation one then has to cover many scales (the box 
size L, \i ~ 0.2 GeV — fOGeV and the lattice cut- 
off a -1 ) simultaneously, the task to reliably match 
the low energy regime with the high energy one, 
i.e. the applicability domain of perturbation theory, 
gets quite complicated. In the present context two 
implementations of such schemes are available, the 
regularization independent approach [Q and the QCD 
Schrodinger functional (SF) [|[ |9|. Whereas the for- 
mer may suffer from the scale hierarchy problem in 
practice, the basic strategy of the SF approach is to 
recourse to an intermediate finite-volume renormal- 
ization scheme, where one identifies two of the before- 
mentioned scales, \i = 1/L, and takes low energy data 
as input in order to use the non-perturbative renor- 
malization group to scale up to high energies || || . 

A schematic view of a non-perturbative compu- 
tation of short distance parameters on the lattice 
along these lines, here in case of the running QCD 
coupling a(/x), is given in the diagram above; the 
same can also be set up for the running quark masses. 
It is important to note that all relations '— are ac- 
cessible in the continuum limit and in this sense uni- 
versal by construction. 



3 Aqcd and M quark)RG i via the SF 

The Schrodinger functional is the QCD partition 
function with certain Dirichlet boundary conditions 
in time imposed on the quark and gluon fields, for 
which a renormalized coupling constant can be de- 
fined as the response to an infinitesimal variation of 
the boundary conditions || . By help of the so-called 
step scaling function, being a measure for the change 
in the coupling when changing the box size L (and 
thus having the meaning of a discrete /3-function) , 
one is now able in the SF scheme to make contact 
with the high-energy regime of perturbative scaling: 

A = lim WW)-' 1 ''^- 1 /^) 

&o = 11/(4tt) 2 , bi = 102/(4tt) 4 . (3) 

Every step during the non-perturbative evolution to- 
wards the perturbative regime has been extrapolated 
to the continuum limit in the quenched approxima- 
tion |l(J , and upon conversion to the MS scheme this 
results in a value for the A-parameter: 

A^= 238(19) MeV. (4) 

An extension of this investigation to the situation 
with two dynamical quarks is already in progress by 
the ALPHA Collaboration. 

In a very similar way, in terms of the current 
quark mass renormalization factor Zp of cq. (^|) re- 
placing the SF coupling to build up another step scal- 
ing function, the scale and scheme independent renor- 
malization group invariant (RGI) quark masses 

M = lim \(2b g 2 (n)y do/2bo m(v)] 

H— >oo I J 

b = 11/(4tt) 2 , do = 8/(4tt) 2 (5) 

were obtained in the same reference. Both evolutions 
are displayed in Fig. 0, and at the scale y,o (leftmost 
point in Fig. |l|b) the matching between the lattice 
regularization and MS via the SF is completed: 

M 

=m = 1.157(12) , Mo ~ 275 MeV . (6) 
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Figure 1: Non-perturbative scale evolution of asF 
and m SF /M computed from simulations of the SF 
in the quenched approximation. The lines represent 
perturbative predictions involving the 2- and 3-loop 
/3-function (a) and 1/2-, 2/2- and 2/3-loop expres- 
sions for the t- and /3-f unctions, respectively (b). 



For the 0(a) improved theory and a massless 
renormalization scheme as utilized here, these results 
can also be summarized as 
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where m is the bare current quark mass, and the 
flavour independent total renormalization factor Zm, 
non-perturbatively known for a range of bare cou- 
plings g in the quenched approximation [|o| , is com- 
posed of an universal part, M/m, and of m/m — 
Za/Zp depending on the lattice regularization. 



4 The strange quark's mass 

In order to illustrate the non-perturbative quark mass 
renormalization just explained in a concrete numer- 
ical application, we first sketch our strategy for the 



computation of light quark masses 1 14 . Their ratios 
are known from chiral perturbation theory (%PT) [ pT| 
as 

£=0.55(4), 
M d M 



24.4(1.5) 



(8) 



with M = |(M U + M d ) Q. Nevertheless there are 
still questions, which might be answered decisively 



only using Lattice QCD. They concern the applica- 
bility of xPT in general, i.e. in how far the lowest 
orders dominate the full result, and the problem that 
the parameters in the chiral Lagrangian (at a given 
order in the expansion) can not be inferred with great 
precision from experimental data alone. This state- 
ment holds in particular for the overall scale of the 
quark masses, which is only defined once the con- 
nection with the fundamental theory, QCD, is made. 
Since the parameters in the chiral Lagrangian (the so- 
called low energy constants) are independent of the 
quark masses, it is important to realize that these 
problems can be dealt with by working with unphys- 
ical — of course not too large — quark masses, where 
it is essential or at least of significant advantage to ex- 
plore a certain range of quark masses. While a deter- 
mination of some low energy constants based on these 
ideas has been recently tested in |l5| , we focus in the 
following on the computation of the renormalization- 
group invariant mass of the strange quark by com- 
bining %PT with lattice techniques. 

To this end, and in the spirit of the considera- 
tions before, we define a reference quark mass M re f 
implicitly through 



m PS (M ief )rl = (m K r ) 2 = 1.5736. 



(9) 



Here 



m PS (M) is the pseudoscalar meson mass 
as a function of the quark mass for mass- 
degenerate quarks, and ro = 0.5 fm and |(m^ + + 
m Ko)| puro QCD = (495MeV) 2 enter the r.h.s. of 
eq. (§). %PT in full QCD relates M rcf to the other 
light quark masses viz. 
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M s + M , 



(10) 



which has been substantiated also numerically in the 
case of quenched QCD |L4|. The remaining task is 
now to calculate M re f from Lattice QCD. 

As the foregoing discussion holds true in mass 
independent renormalization schemes too, one ar- 
rives by virtue of the PCAC relation applied to the 
vacuum-to-pseudoscalar matrix elements at the cen- 
tral relation 



2r M re f 
R 
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PS 



(11) 



where Zm is the flavour independent renormalization 
factor of the previous section, which directly leads 
to the RGI quark masses, being pure numbers and 
not depending on the scheme. By means of numer- 
ical simulations of the SF in large volumes of size 
(1.5 fm) 3 x 3fm, the ratio R/a and the meson mass 
mpsa can be computed accurately as a function of 
the bare quark mass and the bare coupling by evalu- 
ating suitable correlation functions (l3[ |l4j . With the 
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values for the scale r^/a from (T(|, a mild extrapola- 
tion yields R/a at the point mpgr 2 , = 1.5736. Then 
the quantity 2roM rc f is extrapolated to the contin- 
uum limit. Both fits are shown in Fig. ^. In view of 
the still significant slope in the latter, we discard the 
point furthest away from the continuum in this ex- 
trapolation as a safeguard against higher order lattice 
spacing effects. Moreover, the analysis was repeated 
for M le { in units of the kaon decay constant, which 
amounts to substitute eq. (Ill]) by 
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1.5736. 



(12) 



Here we observe a weaker lattice spacing dependence. 
The final results of these analyses 



2r M Iei = 0.36(1) 

vq— 0.5 fm 



2M, 



ref 



w. - a87(3) 

(-Fk)r=160 MeV 



2M ref = 143(5) MeV 



2M ref = 140(5) MeV 



are completely consistent with each other. But, as 
also pointed out in that reference, the assignment 
of physical units is intrinsically ambiguous in the 
quenched approximation. Consulting e.g. the recent 
results of the CP-PACS Collaboration @, roughly 
10 % larger numbers would be obtained, if the scale 
ro were replaced by one of the masses of the stable 
light hadrons. MS masses for finite renormalization 
scales [i can be obtained through perturbative conver- 
sion factors known up to 4-loop precision. A typical 
result is 

mp(2 GeV) = 97(4) MeV , 



(13) 



where the uncertainty in A^L, eq 



, entering the re- 
lation of the running quark masses in the MS scheme 
to the RGI masses, eq. (||), and the quark mass ra- 
tios from full QCD chiral perturbation theory, eq. (||), 



were taken into account 14 



A compilation of lattice results on the strange 
quark mass in the quenched approximation can be 
found e.g. in |is|| . Most of these differ in the Ward 
identity used and in whether non-perturbative renor- 
malization and a continuum extrapolation has been 
performed or not; also systematic errors often are not 
estimated uniformly either. Our result ( |l3| ) includes 
all errors except quenching. Finally it is interesting 
to note that, as reported by the CP-PACS Collabora- 
tion in their comprehensive study about simulations 
with two dynamical flavours ]Ts| , dynamical quark ef- 
fect appear to decrease the estimates for the strange 
quark mass by ~ 20 % or less. 
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Figure 2: Extrapolations of the ratio R to the kaon 
mass scale and, in units of ro, to the continuum limit. 



5 The static-light axial current 

Let us turn to another example, where a scale and 
scheme dependent renormalization is encountered, 
i.e. the matrix element (0|(-i4r) ai |B(p)) = ip^F-g de- 
scribing leptonic B-decays in the theory with heavy 
quarks. It involves the renormalized axial current, 
(^r)m = ZpJ>l[il5d, and the decay constant Fb, 
which is by its own an interesting quantity for a 
first principles computation on the lattice. Since 
mb — 4 GeV 3> Aqcd implies large discretization er- 
rors of O ((omb) 2 ) , a direct treatment assuming a rel- 
ativistic b quark is difficult on the lattice. Therefore, 
in the first place one may restrict to an effective the- 
ory, one possibility being the static approximation, 
where the b quark is taken to be infinitely heavy. 

As at the end we want to relate the physical ma- 
trix element $ at a scale fj, = mb, 



Aqcd 



(14) 



to the one determined on the lattice at some matching 
scale /io, a crucial ingredient is its (scale and scheme 
independent) renormalization group invariant coun- 
terpart 



RGI 
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6 = H/(4^) 2 ,7o--l/47r 2 (15) 
to be passed into the factorization 

$rgi 
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(16) 
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Figure 3: Non-perturbative running of $/$rgi with 
the energy scale in the static approximation, com- 
puted in the SF scheme, compared to perturbation 
theory based on several combinations of orders, to 
which the /3- and 7-functions have been evaluated. 



As already anticipated in the notation, the further 
strategy is basically analogous to that explained when 
considering the coupling and the quark masses: we 
again adopt the SF framework and invoke an appro- 
priate step scaling function, while everything is meant 
in the static approximation now. 

The definition of the renormalized static ax- 
ial current and the step scaling function, together 
with the 2-loop anomalous dimension, has recently 
been worked out perturbatively [ ^o| . The prelimi- 
nary status of the outcome of the corresponding non- 
perturbative investigation by numerical simulations 
of the SF in the quenched approximation is depicted 
in Fig. ||. The leftmost factor in eq. ( fijil ) is then 
supposed to be inferable in that region, where per- 
turbation theory is feasible again. 



6 Conclusions 

Numerical simulations on the lattice can be applied 
to renormalization problems in QCD. In particular, 
the Schrodinger functional scheme offers a clean and 
flexible approach to deal with the accompanying scale 
differences. As a consequence of good control over 
statistical, discretization and systematic errors, non- 
perturbative coupling and quark mass renormaliza- 
tion can be performed with confidence, and solid re- 
sults for A^L and m^ s with high precision of the 
order of a few % were reached in the quenched ap- 
proximation. Similar ideas are now carried over to 



the heavy quark sector of QCD, where first steps to- 
wards a computation of renormalization group invari- 
ant matrix elements in the static approximation are 
under way. 

The presented concepts will be valuable also for 
full QCD. Despite more powerful (super-)computers 
continuously being developed, a quantitative under- 
standing of dynamical sea quark effects is a great 
challenge which, albeit in sight, still demands for 
much effort on the theoretical as well as on the tech- 
nical/implementational side of Lattice QCD. 
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